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S ' Abstract 

, Context. Among late-type red giants, an interesting change occurs in the structure of the outer atmospheric layers as one moves to later spectral 
' ' ■ types in the Hertzsprung-Russell diagram: a chromosphere is always present, but the coronal emission diminishes and a cool massive wind steps 

^ , Aims. Where most studies have focussed on short-wavelength observations, this article explores the influence of the chromosphere and the 
' wind on long-wavelength photometric measurements. 

' Methods. The observational spectral energy distributions are compared with the theoretical predictions of the MARCS atmosphere models for a 
^SJ sample of nine K- and M-giants. The discrepancies found are explained using basic models for flux emission originating from a chromosphere 
T-H ' or an ionized wind. 

j^-^ , Results. For 7 out of 9 sample stars, a clear flux excess is detected at (sub)millimeter and/or centimeter wavelengths, for the other two only 
^— ^ ■ observational upper limits are obtained. The precise start of the excess depends upon the star under consideration. The flux at wavelengths 

' shorter than ~ 1 mm is most likely dominated by an optically thick chromosphere, where an optically thick ionized wind is the main flux 

, contributor at longer wavelengths. 

IL" ' Conclusions. Although the optical to mid-infrared spectrum of the studied K- and M-giants is well represented by a radiative equilibrium 
. 5^ ' atmospheric model, the presence of a chromosphere and/or ionized stellar wind at higher altitudes dominates the spectrum in the (sub)millimeter 
, and centimeter wavelength ranges. The presence of a flux excess also has implications on the role of these stars as fiducial spectrophotometric 
5_j calibrators in the (sub)millimeter and centimeter wavelength range. 
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1. Introduction 

Several studies in the past three decades (see, e.g. , 
Linskv&Haisch'' 19791; lAvres et al.|[l98ll iHiinsch etal.lll996t 
Haisch et al.„199Q) have revealed the existence of dividing lines 
in the cool half of the Hertzsprung-Russell (HR) diagram, 
where the giants and supergiants reside. These dividing lines 
are based on differences in the physics of the outer atmospheric 
layers. To the blue side of the dividing lines, the late-type stars 
are surrounded by chromospheres and coronae. To the red side 
the stars also possess chromospheres, but in combination with 
a cool stellar wind. At the tirne of the introduction of the di- 
viding lines ( Linskv & Haisch 19791; Avres et al. 1981 ) there 
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were no observational indications for the presence of a corona 
on the red side, nowadays there is evidence for some coronal 
emission, alth ough much weake r than on the blue side of the 
dividing hnes ( lAvres et al ] |l997l) . 

Most studies of these outer stellar layers have focussed on 
X-ray and UV observations. But the far-infrared (FIR) contin- 
uum can also be used as a probe of the layers of the atmosphere 
for solar and cooler stars. As the primary infrared (IR) contin- 
uum opacity, coming from free-free processes, increases with 
the square of wavelength, we see emission from increasingly 
outer lying layers the longer the wavelength we observe. 

In this study, we use (sub)millimeter and centimeter wave- 
length observations to gain further insight into the outer struc- 
ture of nine giants of spectral type K and M. The selected stars 
are 'standard' stars, used in the calibration pedigree of many IR 
spectroscopic and photometric instruments. All of them belong 
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Table 1. Some characteri stics of the sample s t ars. The values 
for Mhoi ar e taken from Decin et al.l (l2003ah . iKashvap et al 



(I1994 and lEggen & StokesI digToT No uncertainties 



given in these articles on the values of Mtoi for l Aur and a Lib. 





spectral type 


Mboi 




a Boo 


K2 IIIp 


-0.90 ± 0.05 


1.23 


«Cet 


M2 III 


-3.09 ± 0.13 


1.66 


aTau 


K5 III 


-1.72 ± 0.06 


1.54 


f3 And 


MO III 


-3.14 ± 0.11 


1.58 


/3Peg 


M2.5 III 


-3.34 ± 0.11 


1.67 


7Dra 


K5III 


-2.07 ± 0.07 


1.52 


/3UMi 


K4III 


-1.71 ± 0.07 


1.51 


t Aur 


K3 II 


-2.4 


1.55 


a Lib 


M3/M4 III 


-3.6 


1.65 



to the group of stars with no or low coronal activity. With the 
launch of the ESA-satellite Herschel, which covers the full 55 
to 672 /im wavelength range, it is of interest to study whether 
these low-activity stars can also be used as calibrators at these 
far-infrared wavelengths. In Sect.|2]the selection criteria for the 
stars are outlined. Data reduction for the different data-sets is 
presented in Sect. [3] Sect.|4]confronts the spectral energy distri- 
butions (SEDs) of the stars with the hydrostatic MARCS atmo- 
sphere models. In Sect.|5] the discrepancies found in the previ- 
ous section are discussed. The conclusions are given in Sect.|6] 



Table 2. Fluxes at 1.2 mm determined from SIMBA observa- 
tions. Also listed are the ideal aperture used and the rms noise 
on the sky background. The given uncertainty does not take the 
uncertainty on the absolute calibration into account. 



target 


flux 


ideal aperture 


rms 




(mJy) 


(arcsec) 


(mJy) 


a Boo 


105.6 ± 6.9 


45 


8.4 


(3 And 


< 40.1 




13.8 


a Get 


58.3 ± 4.1 


55 


6.0 


/3Peg 


29.2 ± 2.5 


20 


6.6 



Catalog (iSmith et al 



and lCohen et al. | (|2005h. pho tometric data in the Johnson's 11- 
color sys tem (lDucatill2002l) . the COBE DIRBE Point Source 



2004 and the 2MA SS All-Sky Catalog 



of Point Sources (ISkrutskie et al.l l2006h . A summary of the 
available data can be found in Tables lA.ll - IA.5I in the on- 



line appendix. In order to study the outer atmospheric layers, 
(sub)millimeter and centimeter data have been obtained with 
(1.) the SIMBA bolometer array at 1.2 mm at the SEST tele- 
scope, (2.) the MAMBO II bolometer array at 1.2 mm on the 
IRAM telescope, (3.) the SHARC II camera at 350 ^m and 
450 fim on the CSO, and (4.) the VLA at 22 and 43.3 GHz. The 
reduction of each of these newly obtained data-sets is shortly 
discussed in the next paragraphs. 



2. Sample selection 



3.1. SIMBA observations 



The stars in our sample are part of a larger set of stan- 
dard stars used for the spectrop hotometric calibration o f near- 
and mid-infrared instruments jDeciii et al. l2003blcllai 12004 



Decin & Erikssonl2007 : Gordon et al .l2007h . They are selected 



for their low chromospheric and coronal activit y (see, e.g., 
Obrien & Lambert 



Wiedemann et al. 199 



ospl 
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in Sect. |5.4| an |5.3l l. The position of our sample stars in the 
HR-diagram indicates that they have low coronal activity, how- 
ever, it also points to the presence of a cool stellar wind. These 
stars are potential candidates to be selected as fiducial calibra- 
tion sources for the PACS (60 - 210 ^tm) and SPIRE (200 - 
670 /xm) instruments, which are onboard the ESA Herschel- 
satellite. Hence, it is necessary to check if the possible pres- 
ence of a chromosphere, corona or stellar wind causes a flux 
excess in the far-infrar ed, although they are no t visible in the 
near-infrared CO lines ([Wiedemann et al ] |1994 . Some charac- 
teristics of the selected stars are given in Table [T] 

3. Observations and data reduction 

To construct the SEDs of each of the standard stars sev- 
eral photometric data points were gathered from the literature. 
We have used the U BVRIJKLMNH Photoelectric Catalogue 
dMorel & Magnenall [l978l), the R evised AFGL (RAFGL) 
Catalogue dPrice &Murdocklll983h . the IR AS catalogue of 



Point Sources (IRAS PSC), Version 2.0 dBeichman et al 



19881). observa tions in the Geneva Photo metric Sy s tem 4 
dRufeneiill989h . radio continuum data from IWendkeii d 1995b 



a Boo, (3 And, aCet and /3Peg were observed with SIMBA 
(2003 July 13 - 15) at 1.2 mm, using the fast-scanning tech- 
nique. The MOPs{3 software developed by R. Zylka was used 
for the data reduction. In a first reduction step, some fundamen- 
tal operations like despiking, opacity coiTection and sky-noise 
reduction are performed on each scan. Once the scans made 
during different nights are assembled, the position of the source 
is more accurately determined, which can be used for base- 
line definition and for improvement of the sky-noise reduction. 
For the absolute calibration, scans of Uranus were used. The 
model for Uranus is the standard model offered by MOPSI, the 
calibration uncertainties are estimated at 15%. After the data 
reduction, fluxes were determined using aperture photometry. 
For each source, the 'ideal' aperture was determined, being the 
aperture with the highest corresponding signal-to-noise ratio. 
A more vast description o f the data reduction and analysis of 
the SIMBA data is given in lDehaes etal1d2007l) . Table|2]shows 
the determined fluxes together with the ideal aperture used and 
the rms noise on the sky background. 

3.2. MAMBO II observations 

Observations at 1.2 mm with MAMBO II were obtained for 
a Boo, lAut, /3UMi, 7Dra, aTau, P And, aCet, /3Peg and 
a Lib (2003 October-November). The reduction was done with 



' Observers Handbook SIMBA, 2003, edition 1.9, 
|http://www.Is.eso.org/lasiIIa/TeIescopes/SEST/htniiyteIescope-instruments/simba/inde" 
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Table 3. Fluxes at 1.2 mm determined from MAMBO II obser- 
vations, a Tau is not listed here, since the uncertainty on the 
measurement was too large to give any restraints. The observa- 
tions of a Boo and a Cet are discarded in the subsequent anal- 
ysis (see Sect. l3.2l for explanation). 



target flux (mjy) 



Q Boo 
t Aur 
/3UMi 
7Dra 
/3And 
a Cet 
pPeg 
a Lib 



20.8 ± 3.5 
5.5 ± 1.6 
12.2 ± 1.6 
10.1 ± 1.3 

23.5 ± 2.7 

23.6 ± 2.5 
29.5 ± 3.2 
12.1 ± 2.0 



an adjusted version of the MOPSI software called MOPSlcH, us- 
ing the standard scripts provided for the reduction of On-Off 
data. These scripts also include standard reduction steps such 
as baseline fitting, despiking, correlated skynoise filtering, etc. 
Flux calibration is done using a default conversion factor pro- 
vided by MOPSIC, the calibration uncertainties are estimated at 
15%. After these reduction steps, scans of the same source are 
combined to give one result. Table[3]lists the determined fluxes. 
For 4 sources both SEST and IRAM data were available. Both 
measurements coincide for /3Peg and the upper limit deter- 
mined from the SIMBA observations is in agreement with the 
flux measured by MAMBO II for /3 And. For a Boo and a Cet, 
the 2 measurements do not agree within the errors. As the 
MAMBO II observations were performed in service mode and 
the log-files state very unstable weather conditions for a Cet 
and altocumulus clouds right after the observation of a Boo, 
we have more confidence in the results from the SIMBA ob- 
servations (which were performed in visitor mode). For a Boo 
and a Cet the MAMBO II observations are discarded for the 
remainder of the article. Since the MAMBO II observations for 
the other objects are in good agreement with the SIMBA data 
(e.g. /3 Peg), with data from the catalogs (e.g. /? And) and since 
measurements at these long wavelengths are scarce, the data 
are retained for the other objects. 

3.3. CSO observations 

Observations of five giant stars were made at 350 /im and 
450 ^m using the SHARC II camera at Caltech Submillimeter 
Observatory on several nights in 2005 and 2008. Standard 
Lissajous scans were used for the stars and calibrators. The 
weather conditions were favorable: clear skies, low humidity, 
and precipitable water vapor in the range 1 — 2 mm. Occasional 
periods of unstable atmospheric transmission appear to have 
been properly accounted for in the data analysis. Instead of us- 
ing the facility 225 GHz radiometer for atmospheric extinction 
correction, we used (for each observation in Table |4| the tight 
correlation between the observed signals from the calibrators 
(in raw V) and the average full DC voltage of the bolometers 



at the time of the observations to calculate a calibration factor 
which was then applied to the target star. The full DC bolometer 
voltage is responsive to the emission from the atmosphere and 
therefore its transparency. In the analysis, the detectors were 
corrected for their slight nonlinear response. The beam size of 
CSO/SHARC II at 350 and 450 is 8.3 ± 0.3arcsec and 
9.8 ± 0.3 arcsec, respectively, and all of the giant stars are unre- 
solved. The caUbrators are unresolved or only shghtly resolved. 

The absolute flux calibration is based on the IWrigh j (1 1 976h 
mode l for Mars and subseque nt planet observations and analy- 
sis by lGriffin & OrtonI (119931) . From this work, the absolute un- 
certainties in the fluxes of Uranus and Neptune are believed to 
be 5 %. Our submullimetre observations of /3 Peg are calibrated 
directly vs. Uranus and Neptune and are assigned a system- 
atic calibration uncertainty of 10 % in Table |4] For the remain- 
ing sou rces, we used s econdary calibrators havin g fluxes tabu- 



lated b y ISandelll(ll994 . G. Sandell (priv. comm.). 



Jenness et al 



(I2002h . the JCMT/SCUBA flux calibration web site (2005 up- 
date), and our own cross calibration work. Our best estimates 
for the secondary calibrator fluxes are given in the table, and 
the target stars calibrated with respect to them are assigned a 
systematic uncertainty of 15%. In several cases, the statistical 
uncertainties are much smaller than the systematic uncertain- 
ties, so these measurements would benefit from an improved 
knowledge of the fluxes of the secondary calibrators. 

3.4. VLA observations 

The VLA (Very Large Array) measurements were taken in two 
bands: the Q-band (6.9 cm) and the K-band (1.3 cm). In all of 
our observations, we observed in the continuum mode, which 
effectively provides measurements of the total intensity (Stokes 
I) with an equivalent bandwidth of 184 MHzH. The obser- 
vations for a Boo were undertaken on two separate occasions 
- on January 6, 1999, and on January 25, 2004, for /3Peg on 
April 21, 2005. For the 1999 experiment, the VLA was in the 
C configuration, with maximum physical antenna separation of 
^3.4 km. At this time, only about half of the antennas were 
equipped with Q-band receivers and during our experiment 12 
were available for this frequency. The other 15 were tuned to 
K-band for simultaneous observations. The a Boo observations 
were part of a larger program to observe possible sources for 
millimeter wavelength flux calibration that time, and as such 
were hmited to only about an hour in extent. For the 2004 and 
2005 experiments, the VLA was in the B configuration, with 
maximum physical antenna separation of 11.4km. A full 6 
hour observation was dedicated to the star at Q-band. 

Subsequent calibration of the data proceeded in the nor- 
mal fashion for VLA data, in the AIPS reduction package 



( http://www.cv.nrao. edu/aips/|. For all data, the absolute flux 



density scale was set with an observation of 3C286, with as- 



http://www.iram.es/IRAMES/mainWiki/CookbookMopsic 



The VLA receivers actually operate in the two orthogonal circular 
polarizations, with 92 MHz bandwidth in each polarization. Since 
for a Boo and (3 Peg, we expect the two circular polarizations to have 
equal intensity, they are combined into a total intensity polarization 
(Stokes I), for an effective increase of 2 to the bandwidth, yielding 
184 MHz equivalent bandwidth. 
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Table 4. CSO observations at 350 fim and 450 ^m. 



Star 


Wavelength 


Flux 


Statistical 


Systematic 


Observing 


Calibrator 




(pm) 


(mJy) 


Uncertainty 


Uncertainty 


Dates 


Fluxes 








(mJy) 


(mJy) 




(Jy/beam) 


qBoo 


350 


601 


35 




2005 May 10-13 


Arp 220 (10.5) 






538 


53 




2008 Mar 1 


Arp 220 (10.5) 






507 


19 




2008 May 28 


Arp 220 (10.5) 






529 


27 


79 


average 




a Boo 


450 


488 


48 




2008 Mar 2 


Arp 220 (5.4) 






/I /in 


1 1 
i i 




inns \nr n 
zuuo Apr / 


Arp zzu 






442 


11 


66 


average 




aCet 


350 


210 


16 


32 


2008 Sep 22 


Vesta (11.3) 


a Get 


450 


110 


28 


17 


2008 Sep 24 


Vesta (7.1) 


aTau 


350 


530 


20 


80 


2008 Sep 17-18 


Vesta (10.8), Pallas (9.5) 


aTau 


450 


304 


39 


46 


2008 Sep 24 


Vesta (7.2), CRL618 (11.8), 














HLTau (10.4) 


/?Peg 


350 


361 


9 


36 


2008 Sep 22 


Uranus (234), Neptune (92), 














CRL 2688 (49) 


/3Peg 


450 


240 


12 


24 


2008 Sep 23-24 


Uranus (169), CRL 2688 (26.8) 


7Dra 


350 


116 


25 


17 


2008 Sep 22,24 


CRL 2688 (49) 



sumed flux densities of L455 and 2.520 Jy for Q- and K- 
bands, respectively. Uncertainties in this flux density scale are 
~10% at Q-band and 5% at K-band. Observations of the un- 
resolved secondary calibrator J1357H-193 were used to remove 
long timescale (minutes) atmospheric and system fluctuations 
in the data. The derived flux densities of J1357H-193 were 0.668 
and 0.835 Jy at Q- and K-bands in 1999, and 1 . 175 Jy at Q-band 
in 2004 (the level of variation is common with these point-like 
calibration QSOs at radio wavelengths). 

Images were then constructed from the visibilities via stan- 
dard AIPS routines. The images were lightly CLEANed (a few 
lO's of components) to remove the sampling pattern of the ar- 
ray from them. The final total flux density was then calculated 
in five different ways; (l.)hy counting up the flux density in the 
CLEAN components; (2.) by taking the peak flux density in the 
image; (3.)hy counting up the flux density around the central 
location in the image; (4.) by fitting a gaussian to the image, 
and taking the peak of that fit gaussian (we do not actually re- 
solve the star); and (5.) by actually fitting the visibilities them- 
selves to find the flux density of a point source near the image 
center. The final estimated flux density is taken as the median 
of these five estimates. The uncertainty is taken as the aver- 
age of the uncertainty from the image and visibility fits. This 
is only the formal uncertainty, systematic uncertainties must be 
considered in addition to this. These can arise from: inaccu- 
rate flux density scale, bad pointing, bad elevation corrections, 
atmospheric decorrelation, other electronics sources. Of these, 
by far the dominant uncertainty is the flux density scale, as the 
others are accounted for in various ways in the calibration. 

Table |5] shows the resultant flux densities and uncertainties 
(formal only) for the VLA observations. The two observations 
of a Boo at Q-band are consistent with each other, and the K- 
band observation in 1999 is also consistent, given the expected 
spectral index. 



Table 5. Final flux densities from VLA observations of a Boo 
and /3 Peg. The given uncertainty does not take the uncertainty 
on the absolute calibration into account. 



target 


date 


frequency 


wavelength 


flux density 






GHz 


cm 


mJy 


a Boo 


1999-Jan-06 


22.46 


6.9 


1.7± 0.2 


a Boo 


1999-Jan-06 


43.30 


1.3 


3.3 ± 0.4 


a Boo 


2004-Jan-25 


43.30 


1.3 


3.34 ± 0.08 


/3Peg 


2005-Apr-21 


43.30 


1.3 


2.49 ±0.12 



4. Comparison between SED and theoretical 
predictions 

The observational SEDs are compared with the theoreti - 
cal predictions of the SOSMARCS code of IPlez et al] ( Il992h . 
which is a refined ver sion of the original MARCS code of 
Gustafsson et al.l ( Il975h . The synthetic spe ctra were computed 
with TURBOSPECTRUM dPlez et al.lll992h . the improved ver- 
sion of the SPECTRUM program. For an ove rview of the co ntin- 
uum and line opacity lists used, we refer to becinl (l2000h . The 
MARCS model atmosphere code is built on the assumptions of 
local thermodynamic equilibrium, spherical or plane -parallel 
stratification in homogeneous stationary layers and hydrostatic 
equilibrium. 

The geometry of the radiation transfer problem for the K- 
and M-giants in our sample was given by spherically symmet- 
ric layers. Since the MARCS model atmosphere only extends 
up to 200 /im, the far-infrared continuum spectrum was com- 
puted by extrapolation from the continuum theoretical spec- 
trum between 50 and 200 /im. We therefore have determined 
the temperature of the flux forming region where = 1, 
with A ranging from 50 to 200 /im. With H~ free-free being 
the main continuum opacity source, subsequent outer cooler 
layers are sampled for longer wavelengths. Using a logarith- 
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Alpha Boo 




Wave 



Figure 1. Temperature of the atmospheric model layer where 
r A = 1 for wavelengths between 20 and 200 fim (black thick 
line) for a Boo. The temperature of the characteristic layer 
where most of the photospheric flux is formed {T(^rx=i)) for 
the full 200 to 7500 yum wavelength range is derived by extrap- 
olation from the 50 to 200 /xm wavelength range (grey line). 

mic extrapolation, the temperature for the characteristic layer 
where most of the photospheric flux is formed {T^rx=i)) for 
the full 200 to 7500 /im wavelength range is determined (see 
Fig. [T]). The continuum flux at each far-infrared wavelength 
point is then approximated by the blackbody flux at the charac- 
teristic temperature B\ (T) scaled with the appropriate angular 
diameter. Since we are in the Rayleigh- Jeans part of the spec- 
trum, the flux value is quite insensitive to the temperature, i.e. 
dBx(T)/dT is small. 



The accuracy and resolution of today's FIR instruments 
remain currently too poor to constrain the importance of 
line veiling in t h e (su b)millimeter range. The study by 
Decin & ErikssonI (2007) and B. Plez in case of the 40 - 
665 ^m spectrum for a Tau (priv. comm.) shows that molec- 
ular line absorption at a resolution of ~1500 is typically 
less than 1 % beyond 150 /xm. We therefore will compare the 
(sub)millimeter observational data with continuum flux predic- 
tions. 

The input pararn eters for the MARCS code were taken from 
Decin et al.l ( 12003 al) unless indicated otherwise in Table |6] In 
the same article a discussion about the uncertainties on these 
parameters can be found. The models were reddened according 
to the value of th e inters tellar extinction derived fro m the model 
of Arenou et a l.' ('1992 1 using the distances from IPecin et aP 



(l2003ah ( 



Ochsenbein & Halbwachs ( 1999). The values of the 



interstellar extinction and the distances are listed in Table |6] 

Fig.|2] shows the photometric data in comparison with the 
MARCS models. For all targets, the theoretical predictions un- 
derestimate the observations in the millimeter and/or centime- 
ter wavelength area. Where an excess is detectable at 1.2 mm, 
the model underestimates the observations by an average of 
25 %. At centimeter wavelengths, the discrepancy amounts to 



an average of 90 %. In the following section, different causes 
for this excess are explored. 

5. Discussion 

5. 1. Proof for a significant flux excess at 1.2 mm 

Fig. |2] shows clear indications for a flux excess at millimeter 
and centimeter wavelengths. To prove the flux excess, both the 
observational and theoretical uncertainties in the atmosphere 
models should first be investigated. 

Observational uncertainties The uncertainties on the obser- 
vations in the millimeter/centimeter wavelength region are typ- 
ically of the order of 15 %. The IRAS-PSC error bars given 
in the catalogue are the statistical la uncertainty values; re- 
alistic absolute calibration uncertainties are lacking for the 
PSC, but are estimated to be 20% {D. Kester, priv. comm.). 
This higher uncertainty was already clear from a comparison 
between the Infrared Space Observatory - Short Wavelength 
Spectrome ter (ISO-SWS) data and t he IRAS-PSC and IRAS- 
LRS data dVan Malderen et al.ll2004l) . Therefore we have used 
an error bar of 20 % on the IRAS-PSC data in our analysis. The 
IRAS-PSC fluxes are also colour corrected. 
The uncertainties on the near- and mid infrared photometry 
were taken from the catalogs mentioned in Sect.[3] 



Theor etical uncertainties As described in lPecin & ErikssonI 
(I2OO7I) . the uncertainty on the FIR continuum flux predictions 
mainly arise from uncertainties on (1.) the estimated stellar 
temperature and (2.) the neglect of some physical processes. 

(1.) In the FIR, the dominant continuous opacity arises 
from free-free absorption, whose absorption coefficients 
are nowadays known at an accuracy of about 1 % for wave- 
lengths beyond 0.5 j im. over the temperatur e range between 
1 000 and lOOOOK dPecin & Erikssonll2007l) . An uncertainty 
in the estimated stellar temperature may give rise to an un- 
certainty on the continuum predictions of up to 4 % for A-M 
giants. 

(2.) Since we are tracing regions high up in the atmo- 
sphere, density inhomogeneities and patchy temperature struc- 
tures may occur. This kind of 3-dimensional structures are 
not dealt with in the 1 -dimensional MARCS model atmosphere 
code. Luckily, the wavelength regions of interest are in the FIR, 
where the sensitivity of the Planck function to the temperature 
is small. Another important physical process not included in 
the MARCS atmosphere code is the presence of circumstellar 
dust and/or a chromosphere or ionized wind. While the latter is 
the topic of this study, the first excess can be excluded from the 
det ailed analysi s of the ISO- SWS data for 7 targets in o ur sam- 
ple |Decin[eLal] llooli) and I Van Malderen et all (|2004 . l Aur 
and cr Lib were not observed by ISO, but good-quality IRAS- 
LRS data exist for both objects. The IRAS-LRS data show no 
sign of flux excess due to circumstellar dust. 

A remark concerning the angular diameters that were used 
to compute the fundamental parameters for the MARCS mod- 
els is in place here. The angular diameters are computed from 
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Table 6. Input parameters for the MARCS code from Decin et al] j2003a ) unless indicated otherwise: the effective temperature 
Toff in K, the gravity log g in cm/s^, the microturbulent velocity £,t in kms~^, the metallicity [Fe/H], the abundances of carbon, 
nitrogen, and oxygen, the ^^C/^^C-ratio and the photospheric stellar angular diameter 6*^ in milliarcseconds. The calculation 
of the angular diameter is discussed in Sect . 15. 11 The table also contains the distances (in pc) and the values of the interstellar 
extinction Av as derived from the model of lArenou et alj ( 119921) . Values, for which no literature values have been found, have 
been assumed on the basis of analogue objects, and are Usted in italics. 







L Aur 


UMi 


TDra 




Sp. Type 


K2 IIIp 


K3 II 


K4 III 


K5 III 


K5 III 




4320 ± 140 


4160 ± 130 


4085 ± 140 


3960 ±140 


o o r" A 1 1 /I A 

3850 ± 140 


log .9 


1.5 ± 0.15 


1.74 ± 0.36 


1 1 r\ /~io 

1.6 ± 0.02 


1.30 ± 0.25 


1.50 ± 0.15 




1.7 ± 0.5 


3.00 ± 0.5 


2 ± 0.5 


2.0 ± 0.5 


1.7 ± 0.5 


|Fe/H] 


—0.50 ± 0.20 


— O.il ± 0.22 


—0.15 ± 0.2 


A C\f\ 1 A OA 

0.00 ± 0.20 


A 1 cr 1 A OA 

—0.15 ± 0.20 


e(C) 


7.96 ± 0.20 


8.35 


8.25 ± 0.2 


o 1 cr 1 A o 

8.15 ± 0.25 


O O C 1 A OA 

8.35 ± 0.20 


e(N) 


7.61 ± 0.25 


8.35 


8.15 ± 0.25 


O Oi^ 1 A OC; 


O O C 1 A O C 

8.35 ± 0.25 


e(0) 


8.68 ± 0.20 


8.93 


8.83 ± 0.2 


8.93 ± 0.20 


8.93 ± 0.20 


12^/13^ 


7 ± 2 


10 


9 ± 2 


10 ± 2 


10 ± 2 




20.74 ± 0.10 


7.05 ± 0.03 


9.03 ± 0.42 


9.94 ± 0.05 


20.89 ± 0.10 


distance 


11.26 ± 0.09 


166.56 ± 33.31' 


39.87 ± 7.97' 


45.25 ± 0.94 


19.96 ± 0.38 




0.73 ± 0.27 


3.6 


2.49 ± 0.92 


1.72 ± 1.02 


2.30 ± 0.85 


Av 


0.01 ± 0.15 


0.00 ± 0.15 


0.00 ± 0.5 


0.03 ± 0.15 


0.03 ± 0.15 




(3 And 


Of Get 


/3Peg 


a Lib 




Sp. Type 


MO III 


M2 III 


M2.5 III 


M3/M4 III 




Tcff 


3880 ± 140 


3740 ± 140 


3600 ± 300 


3634 ± 110- 




logs 


0.95 ± 0.25 


0.95 ± 0.25 


0.65 ± 0.40 


0.9 ± 0.31 - 




Ct 


2.0 ± 0.5 


2.3 ± 0.5 


2.0 ± 0.3 


3.1 ± 0.5^ 




[Fe/H] 


0.00 ± 0.30 


0.00 ± 0.30 


0.00 


0.00 




e(C) 


8.12 ± 0.30 


8.20 ± 0.30 


8.20 ± 0.40 


8.23 ± 0.04^ 




e(N) 


8.37 ± 0.40 


8.26 ± 0.40 


8.18 ± 0.40 


8.15 ± 0.05^ 




e(0) 


9.08 ± 0.30 


8.93 ± 0.30 


8.93 ± 0.40 


8.93 





g_^2 10 ±2 5±3 10 

13.03 ± 0.06 12.34 ± 0.06 16.43 ± 0.08 11.00 ± 0.05 

distance 61.12 ± 2.84 67.48 ± 3.78 61.08 ± 2.69 90.80 ± 18.16' 

Mg 2.49 ± 1.48 2.69 ± 1.61 i-^^+f^^l 1.5 

Av 0.06 ± 0.15 0.06 ± 0.16 0.03 ± 0.15 0.20 ± 0.17 



MMcWilliaml ( ll99(]|) . ^ Ijudge & Stencell l ll99lh . ^ lTsuiil ( ll99lh . '^ lAoki & Tsuiil ( ll997h . '^ lOchsenbein & Halbwachsl ( Il999h 



Selby or TCS K-band photometry. For f3 Umi, we have used 
the Johnson K-band magnitude of —1.22 (Faucherre et al. 
1983), yielding a magnitude of —1.276 in the Selby system. 
Zeropoints are calculated using the Kurucz theoretical spec- 
trum of Vega, taking into account the observed near-IR excess 
of Vega (Absil et al. 2006). For the Selby photometric sys- 
tem we obtain a zeropoint of 4.0517 10"^" W/m^/^m, for TCS 
4.4506 10"^" W/m^/^m. An uncertainty of 0.01 mag in the K- 
band photometric data (0. 1 mag for (3 Umi) is propagated in the 
computation of the uncertainty on the angular diameter. At ev- 
ery wavelength, the observed angular diameter represents the 
apparent diameter of the stellar surface where ta ~ 1. Since 
at the longer wavelengths, we are tracing layers that lie further 
and further outwards, it is expected that the angular diameter 
increases with increasing wavelength. If the angular diameter 
in the millimeter and centimeter area is considerably larger than 
the assumed value, the MARCS models will underestimate the 
flux in this wavelength area. 



To investigate this, we derived the change in height of the 
continuum forming layers with increasing wavelength for the 
MO giant (3 And. The layer where t^oss — 1 (with Tross the 
Rosseland optical depth) defines the stellar radius, being in case 
of P And i?* = 6.12 lO^^ ^m. The flux at 100 ^lm is formed 
at + 6.90 10^" cm, for 150 fim at i?, + 7.89 10^" cm and 
for 200 /um at i?* + 1.70 10^^ cm. From the MARCS model we 
derive that the flux at 7 cm (this is the longest wavelength for 
which we have observations) is formed at i?* + 2.23 10^^ cm, 
which corresponds to an increase in radius of 3.64%. For the 
other sample stars, comparable numbers are found. This in- 
crease in angular diameter is insufficient to explain the ob- 
served excess. 

In general, the uncertainties on the theoretical flux predic- 
tions are in the order of 5 to 10 %, excluding the effects of a 
chromosphere or ionized wind. The observed flux excesses are 
hence not caused by inaccuracies in the modeling, but are due 
to physical processes in the stars. 
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Figure 2. Comparison between the photo- 
metric data (asterisks) and the continuum 
theoretical MARCS spectrum (full line) for 
the nine sample stars. If several observations 
are available at the same wavelength, only 
the maximum and minimum flux value were 
plotted, except at (sub)millimeter and cen- 
timeter wavelengths as this wavelength re- 
gion is of particular interest here. Most of 
the error bars fall within the symbols for the 
data. A reversed triangle represents an upper 



5.2. Brightness temperature 

Fig. [3] provides another window at studying the flux excess. It 
shows the brightness temperature over the full 5 /.im to 7 cm 
wavelength range. The brightness temperature is defined as the 
temperature of a black body that gives the same flux as the 
model atmosphere at the i ndicated wavelength, and can be writ- 
ten as ( ICohen et al ] l2005l) 

14387.75/A 



Tb{X) 



In 1 



733.4090 0i 



(1) 



where is the observed flux in Jy, the angular diame- 
ter in milliarcseconds, A the wavelength in /im, and Tb{X) the 
brightness temperature in K. 

For comparison also the brightness temperatures from the 
theoretical models are plotted in Fig. [3] The uncertainty on 
the data and on the angular diameter (see previous section) 
has been propagated to determine the error bars on the ob- 
servational brightness temperatures. The differences between 
the theoretical brightness temperatures (as derived from the 
MARCS predictions) taking the angular diameter uncertainty 



into account are indiscernible. The uncertainty on the effec- 
tive temperature has the highest influence on the theoretical 
brightness temperature predictions. This is illustrated in case 
of [3 UMi in Figs.[3]and|4l In the approximation for long wave- 
lengths, the formula for the brightness temperature shows that 
Tb{\) oc l/6'l). Since the angular diameters at long wave- 
lengths might be underestimated when an extra component be- 
sides the photosphere is present, the brightness temperatures in 
Fig.[3]should be regarded as upper limits. 

In Fig. [3] one can clearly see the wavelength region where 
the flux excess starts for each of the sample stars. For /3UMi, 
the observation at 60 /im still agrees with the model within 
the error Moreover, also the IRAS-LRS and ISO-SWS data of 



(3 UM i agree with the model predictions jVan Malderen et al 



20041) . The flux at 100 i^im is in excess of the model. The same 
holds for L Aur. Note that the IRAS flux at 100 ^m for l Aur is 
of lesser quality than the other IRAS points used in our anal- 
ysis. But recently a weak IR flu x excess at 70 ^m h as been 
found in the Spitzer-MIPS data dOordon et al.1 12007|) . which 
confirms our finding that the excess starts somewhere between 
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Figures. Brightness temperature Tb in 
function of the wavelength for the 9 sam- 
ple stars between 5 fim and 7 cm. The full 
line indicates the brightness temperatures 
derived from the MARCS model, the aster- 
isks show the brightness temperatures de- 
rived from the observations. A reversed tri- 
angle represents an upper limit. The error 
bars on the observational data take the un- 
certainty on the observations and on the an- 
gular diameter into account. The PACS and 
SPIRE wavelength ranges are indicated by 
a dotted line, to facilitate comparison to 
the wavelength region were a flux excess is 
seen. As an illustration of the influence of 
the effective temperature, three models with 
different T^ft are shown for /3UMi. Fig.|4] 
gives a clearer view of these models. 



60 /im en 100 /iin. For a Boo, the model and the data coincide 
for wavelengths up until 350 /im. The measurement at 450 
lies above the model prediction and the data at longer wave- 
lengths are all clearly in excess of the predictions. For a Cet 
the excess seems to start between 450 fim and 1.2 mm, as the 
average flux at 1.2 mm lies well above the model. For /3And 
it starts between 100 //m and 1.2 mm; unluckily no measure- 
ments are available between 100 /im and 1.2 mm. For aXau 
and P Peg, the excess starts at longer wavelengths, as the fluxes 
until 1.38 mm (respectively 1.2 mm) are still in accordance 
with the predictions. For 7Dra and crLib all available data co- 
incide with the model predictions, including the measurements 
at 1.2 mm. However, there are no observational data at longer 
wavelengths and the upper limits at centimeter wavelengths are 
such that they do not exclude an excess. 

It was already argumented that the increase in angular 
diameter at longer wavelengths does not lead to significant 
changes in the MARCS predictions. However, possible extra 
components of the stellar atmosphere not taken into account by 
the MARCS-code, such as a chromosphere, can be extended in 



volume in comparison to the stellar photosph ere. Proof for this 
kind o f extension can be found in for example lDrake & Linskv 
(Il986 ). They studied observations at 2 and 6 cm of, a.o., a Boo. 
They treated the radiation at radio wavelengths as originating 
from an optically thick ionized wind (see also Sect. |5.5l l, for 
which they calculated the half-power radius (meaning that half 
of the radio emission originates from within this radius). For 
a Boo this half -power radius at 2 cm corresponds to the stellar 
radius, but at 6 cm this radius had increased to ~ 1.7 i?,. In the 
next sections, we will elaborate on the possibility that the flux 
excess arises from the presence of a chromosphere (corona) or 
from thermal emission in an ionized stellar wind. 



5.3. Influence of a chromosphere at NIR 
wavelengths 

In 1994, Wiedemann et al. studied the fundamental vibration- 
rotation lines of CO (at /im) in a set of late-type stars. 
The CO Au = 1 lines are remote sensors for the thermal con- 
ditions in the outermost layers of the atmosphere. In particular, 
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og A [firn] 

Figure 4. As an illustration of the influence of the effective 
temperature, three models with different Tcs are shown for 
P UMi. The effective temperature equals 3885, 4085 or 4285 K. 
A different scale then in Fig.|3] is used to enhance the re- 
gion where the models differ. The observation at 6.14 cm is no 
longer visible here. 



flie strongest CO Av = 1 lines occur at or above the tem- 
perature minimum in ch romospheric solar and stellar models. 
Wiedemann et al.l (119941) found observational evidence in favor 
of a 'thermal bifurcation' model for the atmospheres of their 
sample stars. It consists of two distinct physical phases that co- 
exist at chromospheric altitudes. One component is controlled 
by molecular cooling and is represented by a radiative equilib- 
rium model atmosphere with CO induced temperature depres- 
sion. The second component is chromospheric and features a 
temperature inversion produced by the deposition of mechan- 
ical energy. Any observed spectrum from an atmosphere with 
thermal bifurcation is to be interpreted as a spatial sum over 
the two types of thermal regions with appropriate geometrical 
we ighting factors. 



Wiedemann et al 



(119941) concluded that the observed in- 
frared CO Av = 1 spectrum at 4.6 fim of one group of stars, 
containing a Boo, a Tau and 7 Dra, is well described by homo- 
geneous radiative equilibrium models. The near-IR CO spec- 
tra for this group of so-called 'quiet' stars indicate that the 
cool regions dominate the stellar surface for heights between ^ 
10-1/2 and 10-2 g/cm^ in mass column density, and have large 
filUng factors. These stars are said to have a 'COmosphere'. 
For a second group of stars, the CO Av = 1 spectrum is poorly 
represented by radiative equilibrium models, and is compati- 
ble with a chromosphere covering the stellar surface homoge- 
neously at these altitudes. After investigation of different chro- 
mospheric indicators, it also became clear that the stars in the 
first group show only little chromospheric activity. 

Three stars from our sample, a Boo, a Tau and 7 Dra, be- 
long to the first group of 'quiet stars' i.e. their CO Aw = 1 
spectrum indicates that the COmosphere dominates the ther- 
mal structure at heights between ~ 10-^/^ and IQ-^ g/cm^ 
in mass column density. For these stars, the spectrum at NIR 
wavelengths is not influenced by their chromospheric activity. 



In the following sections, we investigate if this remains true at 
longer wavelengths. 

5.4. Coronal, transition region and wind dividing 
lines 



Linsky &Hai"schl(ll979l) introduced a dividing line (further de- 
noted by DL) in the cool part of the HR diagram on the basis 
of ultra-violet spectra of late-type stars. Stars to the blue side 
of the 'transition region DL' were termed 'solar- type', as they 
showed spectral lines formed at temperatures of 5 x 10"^ — 2 x 
10^ K, indicative of chromospheres, transition regions and by 
implication unseen coronae at even hotter temperatures. Stars 
to the red side are of 'non-solar type': they only exhibited lines 
formed at temperatures below 10 000 — 2 000 K, indicative 



of chromospheres only. lAyres et al.l (119811) attempted to ob 



serve soft X-ray emission from late-type stars, being a sig- 
nature of stellar coronae (T > 10^ K). They found an 'X-ray 
DL' roughly coinciding with the 'transition region DL' from 



I979I): only stars t o the blue side were de- 



Stencel & MuUanI (Il980h studied the mor- 



Linskv & Haisch 
tected in X-rays, 
phology of Mg II h and k resonance lines and they also found a 
similar distinction in position in the HR-diagram between stars 
with a low-velocity wind in their chromosphere (to the red side 
of the line) and stars without (to the blue side). 

Since the first discovery of the DLs, several authors have 
confirmed the dichotomy in late-type giant atmospheres, but 
more sensitive observations have a lso revised the location of 
the different DLs in the HR diagram. Hiinsch et al. d 19961) place 
the 'X-ray D L' at B — V = 1.2 for luminosity class III giants. 



According to iReimers et al. I (ll996l) . the 'wind DL' runs verti 
cally at B - V w 1.45, foiB-V < 1.45 i t runs nearly 
horizontal to i? - V" « 1.0 at Alboi « -2.8. And lHaisch et al 



I990I) found that the latest occurrence of emission lines from 



C IV or Si IV, indicative of a transitio n region, occ urs at K4 III, 
corresponding to B — V = 1.385 (lGraylll992l) . This places 
all of our eight stars with luminosity class III in the category 
of late-type stars with a chromosphere and a cool wind (see 
Table [T] and Fig.|5]l. a Boo, although located to t he blue side 



of the wind DL, is known to posses a cool wind dAvres et al 



,1982) and is regarded as an archetype of a no n-coronal star. 



The K 2 II giant l Aur was detected in X-ra vs bv IReimers et al 
d 19961) and shows C IV and Si IV emission ("Haisch et al."l99i 



but also circumstellar Ca II H+K lines (Reimers et al. 199q). It 
is a so-called 'hybrid' giant. Hybrid giants are located to the red 
(i.e. non-coronal side) of the 'transition region DL', but show 
the existence of both transition r egion plasma and large mass- 
loss rates (10-^ - 10-^° Mo/yr) (lHarpeijll992h . 

It should, however, be noted that the distinction between 
coronal and non-coronal giants is not so clear-cut. In a number 
of recent articles, evidence is given that all giants show some 
level of X-ray emission when observed with sufficient sensitiv- 
ity. On the other hand, all observations confirm the significant 
drop in X-ray emission around early K spectral type for lumi- 
nosity class III. 

X-ray emission from cool stars is linked to the confine- 
ment of hot coronal matter in magnetic structures. An under- 
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Figure 5. The position of the different DLs in the HR diagram, 
together with the eight luminosity class 111 gian ts fr o m our sam- 
ple. The 'X-ray DL' is d erived by [Hi insch et al.l ( 1996h . the 
'transit ion region DU by iHaisch et al. |(ll990j), and the 'wind 
DL' by lReimers et al.l (119961) (see text for more details). 



standing of the origin of these magnetic structures is closely 
related to the physical explanation for the existence of the DLs, 
both of which are still under discussion. Most authors seem to 
agree on the fact that a change in surface magnetic field topol- 
ogy is responsible f or the existence of the DLs. According to 
Rosner et al.l (Il995h . the field topology of a red giant changes 
from a large scale organized and closed configuration binding 
coronal gas, to a largely open magnetic field giving rise to a 
massive cool wind, as the star evolves along the RGB from 
the left of the DLs to the right. The transition in topology is 
ascribed to a change in the origin of the field: as the stellar 
rotation drops below a critical value, the spin-catalyzed dy- 
namo gives way to a field generation mechanism requiring fluid 
turbulence as found in the convection zone. It has also been 
suggested that the magnetic flux tubes that rise up from un- 
derneath the convection zone to the stellar surface where they 
form large scale coronal loops, become trapped in the convec- 
tive enve lope as the convection zone deepens to the right of 
the DLs ([Holzwarth & Schiissler 2001). However, Ayres et al 



d 19971120031) have found evidence in "non-coronal" giants that 
coronal loops can still rise to the stellar surface: the loops ex- 
tend beyond the cold molecular layer just above the stellar 
photosphere, but are at least partially submerged in the chro- 
mosphere/COmo sphere, where the coronal X-rays are attenu- 
ated by overlying material. It is also still unclear if stars evolv- 
ing along the RGB cross the DLs (as was a.o. postulated by 
Rosner et al.l 1995h). or if the evolution tracks run parallel to 



them ( Hiinsch 



99m 
et al.1 



1996) 



In this last scenario, the difference 
in X-ray activity on either side of the DLs would be due to a 
different rotational history of each star implying a difference in 
spin-catalyzed dynamo strength. 

In the following sections we will investigate if the observed 
flux excess in the (sub)millimeter and centimeter wavelength 
range in our selected late-type giants can be explained by the 
presence of a chromosphere (corona) and/or by thermal emis- 



sion from an ionized stellar wind. As we shall show in Sect.| 
and Sect. l5.6l the radiation coming from a chromosphere or 
from an ionized wind exhibits a different frequency depen- 
dence cc I/". The spectral indices a as determined from 
the available observations will be compared to the theoretical 
predictions to determine the cause of the flux excess. 

5.5. Thermal emission from an ionized wind 

Stars with an ionized wind emit an excess of continuum emis- 
sion at long wavelengths, i.e. from the IR to the radio region, in 
addition to the black body flux emission. This excess flux is due 
to free-free emission or Bremsstrahlung from the wind. To de- 
rive the wavelength dependence of the the rmal emission from 
an ionized wind, we will use a model from lOlnoiil d 1 975b . This 
article gives analytic expressions for the flux originating from 
a stellar wind, assuming a homogeneous, spherical geometry 
with a uniform electron temperature and with H II as the only 
constituent. In reality, the hydrogen in the winds of these late- 
type stars will be only partially ionized (IDrake et al .Il987h . The 
free-free absorption coefficient per unit m ass k^/ in cm^ g~^ of 
an ion ized gas at long wavelengths is (e.g. lLamers & Cassinelli 
1999h 



P 

where is the square of the charge of the atoms, rig and rii 
are the electron and ion densities in cm~^, p is the density in 
g cm~^ and is the gaunt factor which will be approximated 
by a power law 

.9, ~1.37r° i35A" o«4 (3) 

where A is expressed in cm. It is clear from these expressions 
that the wavelength dependence of the emitted flux will be the 
same, regardless of the main contributor to the flux, be it H~ff 
or H I ff . 



Model V from lOlnonI ( 119751) . the truncated power law dis- 
tribution is of particular intrest here. It assumes an electron 
density distribution with a homogeneous sphere in its centre: 
Tie oc r^^ for r > R, for r < R is constant. Using Eq. (|2|i 
in his expressions, the model predicts cx j/O.eii jj^g opti- 
cally thick limit and Fi, cc jy^o "84 optically thin limit. 
These approximations can be generalised to density distribu- 
tions He OC with j3 > 1.5. A value for /3 differing from 
2 can be caused by a non-constant velocity distribution in the 
wind. This is a very plausible explanation if the radio emission 
originates from the wind acceleration zone. A decreasing (in- 
creasing) fractional ionisation rate (this is the number of free 
electrons per neutral hydrogen atom) with radius can also lead 
to a higher (lower) value for (3. An extreme case of this sce- 
nario would be the existence of an outer cutoff radius ro to 
the ionized portion of the wind. Although doubtful for the tar- 
gets in our study, this cutoff might be caused for example by 
the formation of dus t at this location in the wind, 'quenching' 
the ionized material (IDrake & Linskvll986l) . The spectral index 
would change for v < i/q, where vq is determined by the cutoff 
radius. For K to mid M cool wind giants, most estimates of the 
dust-formation region (if present at all) indicate ro/i?* ~ 10. 
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It can be shown that the s pectrum is only influen ced by this 
transition for A > 30 cm (IDrake & Linsky||l986l) . No obser- 
vations beyond this wavelength are being used in this article, 
hence such a spatial restriction of the ionized region is of no 
importance for our discussion. 

In the optically thin case the wavelength dependence is not 
influenced by the value for /3. In the optically thick case we 
have 



Fi, oc v"' with a 



2/3-3.084 
/3-0.5 



(4) 



For (3 = 1.5 we find a = —0.084, which is the same frequency 
dependence as in the optically thin case, a goes asymptotically 
to 2 as /3 ^ +0O, but we do expect (3 to fluctuate around 2. 

5.6. Chromospheric emission 

The continuum radiation from a chromosphere will be mainly 
free-free emission from H~ and H I. The flux can be written as 



Fx = K{fpBx{T)V 



(5) 



in the optically thin case, where V is the volume of the 
emitting region and B\ is the Planc k function. Th is expres- 
sion can be derived f rom E g. (1) in lOlnonI (Il975h . see also 
Skinner & Whitmord (Il987h . With the use of Eq. © and the 
Rayleigh- Jeans approximation for the Planck function where 
B\{T) oc A^**, we find a wavelength dependence of F\ oc 
^-1.916 ^^gg optically thick chromosphere, we are 

looking at a black body with a temperature equal to the elec- 
tron temperature at the layer where t\ = 1. In that case we find 
Fx oc A~*. In the above approximations, the chromosphere is 
treated as a homogeneous region, with uniform densities and 
electron temperature. 

The model above is also applicable to thermal emission 
from a corona. This implies that the same wavelength depen- 
dence will be found for corona and chromosphere and that 
our analysis will not allow to discriminate between these two 
sources of free-free emission. 

5. 7. Spectral indices for our program stars 

Fig.|6] shows the true flux excess (the theoretically calcu- 
lated flux is already subtracted from the data) at millime- 
ter/centimeter wavelengths. We have chosen to plot X^F\ in 
function of wavelength on a logarithmic scale, since the data 
then follow a horizontal line in case of an optically thick chro- 
mosphere. To determine the spectral indices, we searched for 
the best fitting line through our data using a least-squares 
method. The starting point for these lines coincides with 
the first wavelength where a flux excess is noticeable (see 
Sect. l5.2l i. except for lAut (see later in this section). For the 
least-squares fit we did not take any upper limits into account. 
Table|2] summarizes which spectral indices are expected for a 
chromosphere and for an ionized wind based on the simplified 
analytic expressions derived above. 

The least-squares fits for aCet and (3 Peg have slopes of 
1.2. This value lies closest to the spectral index expected for an 





optically thin 


optically thick 


chromosphere 
ionized wind 


A*Fa (X A^ i 


X^Fx cc A" 



Table 7. Summary of the wavelength dependencies of the 
flux derived for simplified models for the chromosphere and 
the ionized wind. The table gives the wavelength dependence 
for an optically thick ionized wind with a density distribution 
Tie OC r^^ with (3 = 2. For (3 = 1.5 we have X'^Fx oc A^'^ and 
for (3 ^ oowe have X'^Fx oc A°. 



optically thick ionized wind, with f3 equal to 2.2. For f3 And, 
a Tau, (3 UMi and a Boo, the slopes have values of respectively 
0.9, 0.8, 0.7 and 0.7. These values lie somewhere in between 
the slope of an optically thick ionized wind (with /3 equal to 
respectively 2.8, 3.1, 3.5 and 3.5) and an optically thick chro- 
mosphere. For t Aur only two fluxes are available in excess of 
the model, among which the IRAS flux at 100 fim which is of 
lesser quaUty. A least-squares fit to these two points leads to a 
slope of —0.5. For a better determination of the slope, we can 
include the IRAS flux at 60 fim, which is of normal quality and 
just barely agrees with the model within the error This fit gives 
a slope of 0.1 (see Fig.|6]l. Hence, only for t Aur, an optically 
thick chromosphere is found. 

For all stars except l Aur and a Boo, the least-squares fits 
are made to data at wavelengths longer than 1.2 mm, because 
only from this wavelength onwards a flux excess was de- 
tectable. In these cases, an optically thick ionized wind is the 
most likely explanation for the observed excess, although the 
density distributions sometimes show quite large deviations 
from the average rig oc r~^. The only star where an optically 
thick chromosphere is seen, the hybrid giant l Aur, has a line 
fitted only to wavelengths shorter than 1.2 mm, because the flux 
excess was already present in this region and no measurements 
at longer wavelengths were available. We therefore propose 
that at shorter wavelengths ( A < 1 mm) an optically thick chro- 
mosphere is being sampled and at longer wavelengths the con- 
tinuum forming layers lie further outwards in the atmosphere 
and the observations show an optically thick ionized wind. It 
is very well possible that for wavelengths slightly longer than 
'^l mm, the optically thick chromosphere is still visible, as this 
would explain the deviating values for /3 in the optically thick 
ionized winds. 

6. Conclusions 

In 7 out of the 9 K- and M-giants examined, a clear flux excess 
at millimeter and/or centimeter wavelengths was found, for the 
other two targets only observational upper limits are available 
at centimeter wavelengths. The selected stars have low chro- 
mospheric and coronal activity and three of them do belong to 
the group of so-called 'quiet' K- and M-giants, where the near- 
infrared CO Av = 1 lines indicate that the CO-cooled regions, 
as predicted by radiative equilibrium models, dominate over the 
chromosphere at altitudes between ~ 10 "^^^ and 10^^ g/cm^ 
in mass column density. On the basis of this study, it seems 
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Figure 6. These figures show the true flux 
excess for the 7 stars in our sample for 
which a clear flux excess is detected (i.e. 
7 Dra and a Lib are excluded, since only 
upper limits are available). The observations 
(with the continuum already subtracted) are 
represented by asterisks. The full lines show 
the linear least-squares fit to the data. If 
there were different measurements at one 
wavelength, the fluxes were replaced by a 
weighted average. Error bars are shown, but 
they often fall within the symbols for the 
fluxes. A triangle represents an upper limit. 
The slope for each target is indicated in the 
upper corner of each panel. 



that for these stars the homogeneous atmosphere models based 
on radiative equilibrium are able to reproduce the CO spectrum 
around 4.6 /im, but clearly fail to reproduce the flux at millime- 
ter and centimeter wavelengths. At these far-IR wavelengths, 
the presence of a chromosphere and ionized stellar wind cause 
a clear flux excess. 

The observed excess at wavelengths shorter than ~lmm 
is most likely to be attributed to an optically thick chromo- 
sphere, where an optically thick ionized wind is being sampled 
at longer wavelengths. The wavelength region where the ex- 
cess starts depends upon the star under consideration. The most 
extreme cases are / Aur and /?UMi, where the excess starts 
somewhere between 60 /^jm and 100 /im. These findings have 
implications for the roles of these standard stars as fiducial cal- 
ibrators for PACS (wavelengths between 60 and 210 /^/m) and 
SPIRE (between 200 and 670 /im). For a Boo the flux excess 
is already present at the SPIRE (but not at the PACS) wave- 
lengths. For a Cet it might be present at SPIRE wavelengths 
and for (i And the excess might already start at the PACS wave- 
lengths, but a lack of observations in these regions makes it 
impossible to indicate the precise start of the excess, a Tau, 
(i Peg, 7 Dra and a Lib show no flux excess in the PACS and 
SPIRE range, but especially for a Lib, only few observations 
are available in the relevant region, l Aur and /? UMi show a 
clear flux excess from 100 /^im onward, and should not be used 
as a calibrator beyond 60 ^m. 
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Appendix A: Summary of the photometric data 
used in this study 

List of Objects 

L Aur' on page 
crLib' on page 
a Boo' on page 
/? And' on page 
aCet' on pageul 
/3Peg' on page|j 
a Boo' on pageH 
t Aur' on page|2 
/3UMi' on page 2 
7Dra' on page|2| 
aTau' onpage|2| 
/3 And' on pagejH 
aCet' on page|2| 
73 Peg' on pageja 
cr Lib' on page|2| 
aTau' on page Ij 
a Boo' on pageJS 
aCet' on page 
73 Peg' on page 
(3 And' on page 
a Boo' on page 
a Cet' on page 
a Cet' on page 
a Boo' on page 
a Boo' on page 
a Cet' on page 
/3Peg' on page_ 
'/3 And' on page|3| 
7? ' on paged 
a Boo' on pagejj 
■/3Peg' on page|3| 
a Boo' on page 3 
a Boo' on page 3 
/3Peg' on pagebl 
a Boo' onpage|4| 
a Boo' on pagels] 
a Tau' on pagejj 
i Aur' on page 
cr Lib' on page_ 
/JUmi' on page|2 
'/3Umi' on pagemi 
/? And' on page|6| 
'(3 And' on pagemi 
^/?UMi' on pageQ 
t Aur' on pageM 
6 Aur' on page|7| 
/?UMi' on page! 
a Boo' onpage|8| 
a Cet' on page|8| 
/? And' on page|^ 
aTau' on page 
/3Peg' on page_ 
7Dra' on page 1 



(7 Lib' on page|8[ 
a Boo' on page|& 
a Boo' on page|8l 
/3 UMi' on pagej] 
a Boo' on pagej] 
a Tau' on pagel 
7Dra' on page[ 
a Boo' on pagej] 
aTau' on pagel] 
7Dra' on page[ 
a Boo' on pagej] 
L Aur' on page 
L Aur' on page 
a Cet' on page 111 
/3 Peg' on pagein 
/3 And' on page hi 
a Tau' on pagell 1 
(3 UMi' on pagefn 
a Boo' on pag ell if 
L Aur' on page H 
L Aur' on page H 
i Aur' on pagell 1 
a Boo' on pag ell ll 
t Aur' on page ll ll 
7 Dra' on pagell2l 
a Lib' on pag ell2l 
tAur' on page ll2l 
(3 UMi' on page[l2 
a Boo' on pag ell2l 
a Cet' on page ll2l 
(3 And' on pagell 
a Tau' on page 12 
/3 Peg' on page 12 
7 Dra' on page 12 
a Lib' on page 12 
crLib' on page 12 
L Aur' on pagell2 
(3 UMi' on pagefl 
a Boo' on pagejj 
L Aur' on page [3 
(3 UMi' on page 
7 Dra' on pageM 
a Tau' on pageD 
/3 And' on pagep 
a Cet' on page 
/3 Peg' on page 
a Lib' on page 7| 
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Table A.l. Photometric data used in this study for the targets a Boo and t Aur. The hterature references are specified at the end 
of Table |A3] 
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Table A.2. Photometric data used in this study for the targets /3 UMi and 7 Dra. The Uterature references are specified at the end 
of Table |A31 
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± 0.00 




1 


4.01 10"^ 


3.91 10" 


10 


± 4.19 10 


-12 


4 


4.01 10"^ 


9 QQ 1 D"^" 


± 


2.38 10 


-12 


4 


4.23 10~^ 


7.79 10" 


10 


± 7.1710 


-12 


4 


4.23 10"^ 


fi 1 7 1 n"^" 


± 4.02 10 


-12 


4 


4.40 10"^ 


1.15 10" 


9 


-|- 


0.00 




1 


4.40 10"^ 


y.DO lu 


± 0.00 




1 


4.48 10~^ 


1.24 10" 


9 


-|- 


1.33 10" 


11 


4 


4.48 10"^ 


1 ni in^^ 


-|- 


8.04 10" 


12 


4 


5.39 10~^ 


2.78 10" 


9 


-|- 


2.99 10" 


11 


4 


5.39 10"^ 


9 1 n~'' 


-|- 


1.88 10" 


11 


4 


5.49 10~^ 


3.01 10" 


9 


-|- 


2.22 10" 


11 


4 


5.49 10"^ 




-|- 


1.18 10" 


11 


4 


5.50 10~^ 


3.06 10" 


9 


-|- 


0.00 




1 


5.50 10"^ 


9 fiQ 1 
z.oy lu 


-|- 


0.00 




1 


5.81 10~^ 


3.66 10" 


9 


-|- 


3.93 10" 


11 


4 


5.81 10"^ 


^10 1 

O.Lit lU 


-|- 


2.54 10" 


11 


4 


7.00 10~^ 


5.28 10" 


9 


-|- 


0.00 




1 


7.00 10"^ 




-|- 


0.00 




1 


9.00 10~^ 


6.68 10" 


9 


-|- 


0.00 




1 


9.00 10"^ 


D.OD lU 


-|- 


3.00 




1 


1.24 


6.10 10" 


9 


-|- 


1.09 10" 


09 


12 


1.24 


/I s/1 1 r\~^ 


-|- 


9.46 10" 


10 


12 


1.25 


6.43 10" 


9 


-|- 


0.00 




6 


1.25 


0.01 lU 


-|- 


0.00 




6 


1.25 


5.64 10" 


9 


-|- 


2.01 10" 


10 


7 


1.25 


K ztR 1 


-|- 


2.20 10" 


10 


7 


1.66 


5.82 10" 


9 


-|- 


9.86 10" 


10 


12 


1.66 


/I 88 1 n^^ 


-|- 


8.09 10" 


10 


12 


2.16 


3.03 10" 


9 


-|- 


5.68 10" 


10 


12 


2.16 


9 70 1 n^^ 


-|- 


3.97 10" 


10 


12 


2.20 


3.09 10" 


9 


-|- 


0.00 




6 


2.20 


9 07 1 r\~^ 


-|- 


0.00 




6 


2.20 


2.85 10" 


9 


-|- 


3.32 10" 


11 


7 


2.20 


9 C0 1 


-|- 


4.24 10" 


11 


7 


3.50 


9.18 10" 


10 




2.94 10 


-11 


7 


3.40 


1 00 1 n^^ 


-|- 


0.00 




1 


3.50 


8.69 10" 


10 




2.94 10 


-11 


7 


3.40 


1 0*^ 1 


± 


0.00 




6 


4.20 


6.02 10" 


10 


± 


1.66 10 


-10 


2 


3.50 


Q 1 S 1 0"^" 


± 


2.45 10 


-11 


7 


4.90 


2.92 10" 


10 


± 


6.45 10 


-12 


7 


3.50 


Q «0 1 fl"^" 


± 


2.45 10 


-11 


7 


4.90 


2.82 10" 


10 


± 


6.45 10 


-12 


7 


4.20 


fi fiO 1 n^^" 


± 


1.82 10 


-10 


2 


1.10 10^ 


4.37 10" 


11 


± 


8.05 10 


-12 


2 


4.90 


9 84 1 n^^" 

Zi.o^ -LU 


± 


6.9410 


-12 


7 


1.20 10^ 


2.82 10" 


11 


± 


5.65 10 


-12 


3 


4.90 


2.74 10"^" 


± 


6.9410 


-12 


7 


1.98 10^ 


7.77 10" 


12 


± 


1.43 10 


-12 


2 


5.00 


3.23 10"^" 


± 


0.00 




6 


2.50 10^ 


3.26 10" 


12 


± 


6.52 10 


-13 


3 


1.02 10^ 


3.81 10"" 


± 


0.00 




1 


2.7410^ 


4.1710" 


12 


± 


1.15 10 


-12 


2 


1.02 10^ 


3.88 10"" 


± 


0.00 




1 


6.00 10^ 


2.1710" 


la 


± 


4.35 10 


-14 


3 


1.1010^ 


3.63 10"" 


± 


6.70 10 


-12 


2 


1.00 10^ 


5.73 10" 


14 


± 


1.15 10 


-14 


3 


1.2010^ 


2.72 10"" 


± 


5.45 10 


-12 


3 


1.20 10^ 


4.00 10" 


17 


± 


1.00 10 


-17 


5 


1.9810^ 


1.62 10"" 


± 


2.99 10 


-12 


2 


1.20 10^ 


3.05 10" 


17 


± 


6.00 10 


-18 


10 


2.5010^ 


3.24 10"^^ 


± 


6.48 10 


-13 


3 


6.1410^ 


5.3710" 


21 


± 


0.00 




5 


6.00 10^ 


2.25 10"" 


± 


4.51 10 


-14 


3 
















1.0010^ 


4.89 10"" 


± 


9.7710 


-15 


3 
















3.5010^ 


9.94 10"" 


± 


2.5710 


-16 


3 
















1.2010^ 


2.52 lO"^'^ 


± 


5.00 10 


-18 


10 
















1.20 10''' 


< 3.0010 


-17 






5 
















2.01 10'' 


< 6.42 10 


-20 






5 
















3.55 10"* 


< 1.2710 


-20 






5 
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Table A.3. Photometric data used in this study for the targets a Tau and (3 And. The literature references are specified at the end 
of Table |A3] 



a Tau 


f3 And 


A 








A 




ref. 


X 




XF 


A 




ref. 


[//nil 




[Wm 








LP J 


[Wm 








3.46 10^^ 


1.85 10" 


10 


± 3.5110 


-12 


4 


3.46 10"^ 


o.oo lU 


± 6.6110 


-13 


4 


3.60 10"^ 


2.93 10" 


10 


± 0.00 




1 


3.60 10"^ 


Q 1 7 1 n^^-'^ 
y . 1 / lu 


± 0.0010 


fOO 


1 


4.01 10~^ 


9.98 10" 


10 




1.90 10 


-11 


4 


4.01 10"^ 


0. 1 1 lU 


± 3.8310 


-12 


4 


4.23 10~^ 


2.09 10" 


9 




3.47 10" 


11 


4 


4.23 10"^ 


fi fin 1 n"^" 

u.uu lU 


-|- 


5.98 10 


-12 


4 


4.40 10"^ 


3.31 10" 


9 


-l- 


0.00 




1 


4.40 10"^ 


1 ns 1 n^^ 

1 .UO lU 


-|- 


0.00 




1 


4.48 10~^ 


3.37 10" 


9 


-l- 


fi 4n 1 n~ 


11 


4 


4.48 10"^ 


1 07 1 n^^ 

L.U t lU 


-|- 


1.31 10" 


11 


4 


5.39 10~^ 


8.08 10" 


9 


-l- 




10 


4 


5.39 10"^ 


9 79 1 n~^ 

^. 1 Z -LU 


-|- 


3.34 10" 


11 


4 


5.49 10~^ 


8.79 10" 


9 


-l- 


1.21 10" 


10 


4 


5.49 10"^ 


9 QA 1 r\~^ 

Z.y^ iU 


-|- 


1.08 10" 


11 


4 


5.50 10~^ 


9.41 10" 


9 


-l- 


n nn 




1 


5.50 10"^ 


^1/11 n~^ 


-|- 


0.00 




1 


5.81 10~^ 


1.07 10" 


8 


-l- 


2.04 10" 


10 


4 


5.81 10"^ 


o.oy lu 


-|- 


4.41 10" 


11 


4 


7.00 10~^ 


1.81 10" 


8 


-l- 


n nn 




1 


7.00 10"^ 


fi 1 9 1 r\~^ 

O. iZ lU 


-|- 


0.00 




1 


9.00 10~^ 


2.71 10" 


8 


-l- 


n nn 




1 


9.00 10"^ 


Q fif^ 1 n^^ 
y.DO lu 


-|- 


0.00 




1 


1.24 


2.66 10" 


8 


-l- 


4 7^ 1 n~ 


9 


12 


1.24 


Q '^'^ 1 n^^ 


-|- 


1.77 10" 


09 


12 


1.25 


2.31 10" 


8 


-l- 


n nn 

u.uu 




1 


1.25 


ft fiR 1 n^^ 


-|- 


2.92 10" 


10 


7 


1.25 


2.38 10" 


8 


-l- 


0.00 




6 


1.25 


q on 1 
y.zu lu 


-|- 


0.00 




6 


1.66 


2.43 10" 


8 


-l- 


Q sn 1 n" 

O.OU -LU 


9 


12 


1.66 


8 SO 1 n^^ 

O.OU lU 


-|- 


1.20 10" 


09 


12 


2.16 


1.53 10" 


8 


-l- 


1 Q7 1 n~ 

i-.U 1 -LU 


9 


12 


2.16 


^ nfi 1 n~^ 


-|- 


7.46 10" 


10 


12 


2.20 


1.14 10" 


8 


-l- 


n nn 

u.uu 




1 


2.20 




-|- 


5.64 10" 


11 


7 


2.20 


1.19 10" 


8 


-l- 


0.00 




6 


2.20 


q no 1 
o.yz lu 


-|- 


0.00 




6 


3.40 


4.33 10" 


9 


-l- 


0.00 




6 


3.40 


1 87 1 n~^ 

1 .o / lU 


-|- 


0.00 




1 


3.50 


3.79 10" 


9 


-l- 


Q 97 1 n" 

iy.j^ 1 -LU 


11 


7 


3.40 


17/11 n~^ 


-|- 


0.00 




6 


3.50 


3.65 10" 


9 


-l- 


9.27 10" 


11 


7 


3.50 


1 1 (-1-9 
1 .UU lU 


-|- 


3.78 10" 


11 


7 


4.20 


2.88 10" 


9 


-l- 


7.96 10" 


10 


2 


3.50 


14-^1 n~'' 

L.^O lU 


-|- 


3.78 10" 


11 


7 


4.90 


1.19 10" 


9 


-l- 


1 "^n 1 n" 

± .OU -LU 


11 


7 


4.20 


Q ^A^ n^^" 

y.u4 lu 


-|- 


2.64 10 


-10 


2 


4.90 


1.17 10" 


9 


-l- 


1 "^n 1 n" 

± .OU -LU 


11 


7 


4.90 


4 7Q 1 n^^" 

4:. t O lU 


-|- 


6.23 10 


-12 


7 


5.00 


1.24 10" 


9 


-l- 


n nn 

U.uu 




1 


4.90 


/I fif^ 1 n^^*^ 

4. DO lU 


-|- 


6.23 10 


-12 


7 


5.00 


1.40 10" 


9 




n nn 

u.uu 




6 


5.00 


/I 7f^ 1 n^^'^ 


-|- 


0.00 




1 


1.02 10^ 


1.52 10" 


10 




n nn 

u.uu 




1 


5.00 


Pi f^i 1 n^^'^ 

0.01 lU 


-|- 


0.00 




6 


1.02 10^ 


1.95 10" 


10 




n nn 

u.uu 




6 


1.02 10^ 


fi Q/i 1 n^^^ 
D.y4 lu 


-|- 


0.00 




1 


1.10 10^ 


1.74 10" 


10 




4.81 10 


-11 


2 


1.02 10^ 


ft 1 1 1 

o.ll lU 


± 


0.00 




6 


1.20 10^ 


1.23 10" 


10 


± 


2.46 10 


-11 


3 


1.10 10^ 


1 .oy lu 


± 


1.40 10 


-11 


2 


1.98 10^ 


3.09 10" 


11 


± 


5.70 10 


-12 


2 


1.20 10^ 


O.UO lU 


± 


1.0110 


-11 


3 


2.50 10^ 


1.3110" 


11 


± 


2.61 10 


-12 


3 


1.98 10^ 


1.12 10"^^ 


± 


2.0710 


-12 


2 


2.74 10^ 


9.55 10" 


12 


± 


2.6410 


-12 


8 


2.50 10^ 


5.84 10"^^ 


± 


1.1710 


-12 


3 


6.00 10^ 


9.79 10" 


13 


± 


1.96 10 


-13 


3 


2.74 10^ 


4.5710"^^ 


± 


1.26 10 


-12 


2 


1.00 10^ 


1.76 10" 


13 


± 


3.52 10 


-14 


3 


6.0010^ 


4.08 10"" 


± 


8.16 10 


-14 


3 


3.50 10^ 


4.5410" 


15 


± 


7.02 10 


-16 


13 


1.0010^ 


8.39 10"" 


± 


1.68 10 


-14 


3 


4.50 10^ 


2.03 10" 


15 


± 


4.00 10 


-16 


13 


1.20 10''^ 


6.25 10~" 


± 


1.00 10 


-17 


5 


1.20 10^ 


1.2710" 


16 


± 


1.50 10 


-17 


5 


1.20 10''^ 


< 1.0010 


-16 






9 


1.38 10^ 


5.62 10" 


17 


± 


1.23 10 


-17 


8 


1.20 10^ 


5.8710"^'^ 


± 


1.1010 


-17 


10 


2.7710^ 


1.5110" 


17 


± 


1.58 10 


-18 


8 


1.3310" 


< 2.02 10 


-17 






5 


2.00 lO'' 


8.98 10" 


20 


± 


1.50 10 


-20 


5 


2.01 10" 


8.3710"^" 


± 


0.00 




5 


2.01 10" 


1.28 10" 


19 


± 


0.00 




5 


2.80 10" 


< 5.3410 


-19 






5 


2.80 10" 


< 5.3410 


-19 






5 


3.55 10" 


1.52 10"^" 


± 


0.00 




5 


3.55 10" 


2.53 10" 


20 


± 


0.00 




5 


6.1410" 


< 1.0310 


-20 






5 


6.1410" 


< 1.3210 


-20 






5 


6.97 10'^ 


< 4.3010 


-19 






5 


6.9710^ 


< 4.3010 


-19 






5 
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Table A.4. Photometric data used in this study for the targets a Cet and /? Peg. The hterature references are specified at the end 
of Table |A3] 





a 


Cet 






/3Peg 


A 








ref. 


A 






\F 


A 




ref. 












LP J 




[Wm 








3.46 10^^ 


3.4410" 


" ± 5.6710 


-13 


4 


3.46 10"^ 


3.36 10" 


11 


± 


1.80 10 


-li 


4 


3.60 10^^ 


5.68 10" 


" ± 0.00 




1 


3.60 10"^ 


5.95 10" 


11 


± 0.00 




6 


4.01 10"^ 


1.98 10" 


± 3.2710 


-12 


4 


3.60 10"^ 


6.23 10" 


11 


± 0.00 




1 


4.23 10"^ 


4.16 10" 


^° ± 5.70 10 


-12 


4 


4.01 10"^ 


2.16 10" 


10 


± 


1.1610 


-11 


4 


4.40 10"^ 


6.49 10" 


± 0.00 




1 


4.23 10"^ 


4.42 10" 


10 


± 


2.32 10 


-11 


4 


4.48 10"^ 


6.82 10" 


^° ± 1.1310 


-11 


4 


4.40 10"^ 


6.99 10" 


10 


± 0.00 




6 


5.39 10"^ 


1.82 10" 


^ ± 3.0010" 


11 


4 


4.40 10"^ 


6.61 10" 


10 


± 0.0010 


f 00 


1 


5.49 10"^ 


1.98 10" 


^ ± 2.01 10" 


11 


4 


4.48 10"^ 


7.17 10" 


10 


± 3.8410 


-11 


4 


5.50 10~^ 


2.02 10" 


^ ± 0.00 




1 


5.39 10"^ 


2.01 10" 


9 


-|- 


1.07 10" 


10 


4 


5.8110"^ 


2.39 10" 


^ ± 3.9510" 


11 


4 


5.49 10"^ 


2.15 10" 


9 


-|- 


1.11 10" 


10 


4 


7.00 10"^ 


4.35 10" 


^ ± 0.00 




1 


5.50 10"^ 


2.24 10" 


9 


-|- 


3 00 




6 


9.00 10"^ 


7.95 10" 


^ ± 0.00 




1 


5.50 10"^ 


2.08 10" 


9 


-|- 


3.00 




1 


1.24 


7.5110" 


^ ± 9.9010" 


10 


12 


5.81 10"^ 


2.53 10" 


9 


-|- 


1.36 10" 


10 


4 


1.25 


7.1110" 


^ ± 0.00 




1 


7.00 10"^ 


5.53 10" 


9 


-|- 


0.00 




6 


1.25 


7.56 10" 


^ ± 2.4010" 


10 


7 


7.00 10"^ 


5.48 10" 


9 


-|- 


0.00 




1 


1.66 


8.2710" 


^ ± 1.3910" 


9 


12 


9.00 10"^ 


1.17 10" 


8 


-|- 


0.00 




6 


2.20 


3.85 10" 


± 0.00 




1 


9.00 10"^ 


1.15 10" 


8 


-|- 


0.00 




1 


2.20 


4.10 10" 


^ ± 4.6310" 


11 


7 


1.24 


1.45 10" 


8 


-|- 


9 ^8 1 n~ 


09 


12 


3.40 


1.36 10" 


" ± 0.00 




1 


1.25 


1.20 10" 


8 


-|- 


4.31 10" 


10 


7 


3.40 


1.53 10" 


^ ± 0.00 




1 


1.25 


1.10 10" 


8 


-|- 


4.31 10" 


10 


7 


3.50 


1.33 10" 


^ ± 4.3610" 


11 


7 


1.66 


1.34 10" 


8 


-|- 


2.22 10" 


09 


12 


3.50 


1.26 10" 


± 4.3610" 


11 


7 


2.16 


8.27 10" 


9 


-|- 


1.17 10" 


09 


12 


4.20 


9.5410" 


± 1.7610 


-10 


2 


2.20 


6.51 10" 


9 


-|- 


0.00 




1 


4.90 


4.08 10" 


± 9.49 10 


-12 


7 


2.20 


6.98 10" 


9 


-|- 


9.16 10" 


11 


7 


4.90 


3.86 10" 


± 9.49 10 


-12 


7 


3.40 


2.42 10" 


9 


-|- 


0.00 




1 


5.00 


3.6410" 


± 0.00 




1 


3.40 


2.56 10" 


9 


-|- 


0.00 




1 


5.00 


4.58 10" 


± 0.00 




6 


3.50 


2.25 10" 


9 


-|- 


5.76 10" 


11 


7 


1.02 10^ 


5.1710" 


" ± 0.00 




1 


3.50 


2.09 10" 


9 


-|- 


5.76 10" 


11 


7 


1.02 10^ 


6.10 10" 


± 0.00 




6 


4.20 


1.51 10" 


9 


± 


4.18 10" 


10 


2 


1.10 10^ 


5.25 10" 


± 1.45 10 


-11 


2 


4.90 


7.20 10" 


10 


± 


1.22 10 


-11 


7 


1.20 10^ 


4.12 10" 


" ± 8.2410 


-12 


3 


4.90 


6.69 10" 


10 


± 


1.22 10 


-11 


7 


1.98 10^ 


8.52 10" 


± 0.00 




2 


5.00 


7.26 10" 


10 


± 


0.00 




1 


2.50 10^ 


4.78 10" 


± 9.5710 


-13 


3 


5.00 


8.26 10" 


10 


± 


0.00 




6 


6.00 10^ 


3.38 10" 


" ± 6.76 10 


-14 


3 


1.02 10^ 


9.94 10" 


11 


± 


0.00 




1 


1.00 10^ 


6.69 10" 


" ± 1.34 10 


-14 


3 


1.02 10^ 


1.24 10" 


10 


± 


0.00 




6 


3.50 10^ 


1.80 10" 


± 3.08 10 


-16 


13 


1.10 10^ 


1.00 10" 


10 


± 


2.7710 


-11 


2 


4.50 10^ 


7.33 10" 


± 2.2010 


-16 


13 


1.20 10^ 


6.82 10" 


11 


± 


1.36 10 


-11 


3 


1.20 10^ 


3.75 10" 


" ± 7.50 10 


-18 


5 


1.98 10^ 


1.95 10" 


11 


± 


0.00 




2 


1.20 10^ 


1.46 10" 


± 2.33 10 


-17 


9 


2.5010^ 


8.29 10" 


12 


± 


1.66 10 


-12 


3 


2.01 10^ 


7.62 10" 


2" ± 0.00 




5 


2.74 10^ 


6.61 10" 


12 


± 


0.00 




2 


3.55 10^ 


6.25 10" 


2" ± 0.00 




5 


6.00 10^ 


5.95 10" 


13 


± 


1.19 10 


-13 


3 


3.56 10^ 


< 1.4310"^** 




5 


1.00 10^ 


1.03 10" 


13 


± 


2.06 10 


-14 


3 


6.1410" 


< 7.8210"^^ 




5 


3.5010^ 


3.09 10" 


15 


± 


3.1710 


-16 


13 


1.31 10'^ 


< 1.0110"^* 




5 


4.50 10^ 


1.60 10" 


15 


± 


1.80 10 


-16 


13 












1.20 10^ 


5.75 10" 


17 


± 


1.25 10 


-17 


5 












1.2010^ 


7.29 10" 


17 


± 


1.22 10 


-17 


9 












1.20 10^ 


7.3710" 


17 


± 


1.3710 


-17 


10 












7.00 10=* 


1.0710" 


18 


± 


5.1410 


-20 


11 












2.01 W 


1.25 10" 


19 


± 


0.00 




5 












3.55 W 


< 3.0410 


-20 






5 












6.14 l(f 


< 9.7710 


-21 






5 
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a 


Lib 








A 










ref. 




Lr^ J 




[Wm"2] 






3 


46 10"' 


1.62 10" 


ii 


± 5.96 10 


-13 


4 


3 


60 10"' 


2.82 10" 


11 


± 0.00 




6 


3 


60 10"' 


2.62 10" 


11 


± 0.00 




1 


4 


01 10"' 


9.69 10" 


11 


± 3.57 10 


-12 


4 


4 


23 10"' 


1.98 10" 


10 


± 6.70 10 


-12 


4 


4 


40 10"' 


3.08 10" 


10 


± 0.00 




1 


4 


40 10"' 


3.28 10" 


10 


± 0.00 




6 


4 


48 10"' 


3.24 10" 


10 


± 1.19 10 


-11 


4 


5 


39 10"' 


9.09 10" 


10 


± 3.35 10 


-11 


4 


5 


49 10"' 


9.77 10" 


10 


± 2.97 10 


-11 


4 


5 


50 10"' 


1.02 10" 


9 


± 0.00 




1 


5 


50 10"' 


1.08 10" 


9 


± 0.00 




6 


5 


81 10"' 


1.16 10" 


9 


± 4.26 10" 


11 


4 


7 


00 10"' 


2.60 10" 


9 


± 0.00 




1 


7 


00 10"' 


2.96 10" 


9 


± 0.00 




6 


9 


00 10"' 


5.35 10" 


9 


± 0.00 




1 


9 


00 10"' 


5.98 10" 


9 


± 0.00 




6 


1 


24 


5.19 10" 


9 


± 7.84 10 


10 


12 


1 


25 


5.65 10" 


9 


± 0.00 




6 


1 


25 


5.20 10" 


9 


± 1.40 10" 


10 


7 


1 


66 


5.70 10" 


9 


± 9.13 10" 


10 


12 


2 


16 


3.35 10" 


9 


± 6.12 10" 


10 


12 


2 


20 


3.20 10" 


9 


± 0.00 




6 


2 


20 


3.06 10" 


9 


± 4.83 10 


11 


7 


3 


40 


1.18 10" 


9 


± 0.00 




1 


3 


50 


1.03 10" 


9 


± 2.60 10" 


11 


7 


3 


50 


9.74 10" 


10 


± 2.60 10 


-11 


7 


4 


20 


6.02 10" 


10 


± 1.11 10 


-10 


2 


4 


90 


3.23 10" 


10 


± 7.76 10 


-12 


7 


4 


90 


3.11 10" 


10 




-12 


7 


5 


00 


3.23 10" 


10 


± 0.00 




1 


5 


00 


3.60 10" 


10 


± 0.00 




6 


1 


02 10' 


3.5110" 


11 


± 0.00 




1 


1 


02 10' 


4.42 10" 


11 


± 0.00 




6 


1 


10 10' 


6.32 10" 


11 


± 1.7510 


-11 


2 


1 


20 10' 


3.55 10" 


11 


± 7.10 10 


-12 


3 


1 


98 10' 


2.1410" 


11 


± 9.86 10 


-12 


2 


2 


50 10' 


3.69 10" 


12 


± 7.39 10 


-13 


3 


6 


00 10' 


3.00 10" 


18 


± 6.00 10 


-14 


3 


1 


00 10^ 


6.8410" 


14 


± 1.3710 


-14 


3 


1 


20 10^ 


3.02 10" 


17 


± 6.75 10 


-18 


10 


3 


5710'' 


< 4.5410 


-19 




5 


6 


14 10* 


< 1.0310 


-20 




5 



' UBVRIJKLMNH Photoelectric Catalogue jMorel & MagnenaJI 19781) , ^ The Revised AFGL (RAF GL) Catalogue jPrice & Murdocklll983l) . 
IRAS catalogue of Point Sources, Version 2.0, * Observations in the Geneva Photometr ic Syste m 4 teufeneilll989h . ^ Radio continuum 



emission from stars ('Wendker"l995'), ^ Stel lar Photometry in Jo hnson's 1 1 -color system l lDucatill2002h 
dSmith et alj2004h . .Cohen et al. ( ,20o3) . 



Dehaes et al 



ll2007h 



COBE DIRBE Point Source Catalog 



2MASS All-Sky Catalog of Point Sources l lSkrutskie et al.ll2006h '^ CSO observations (this ai'ticle) 



'° IRAM observations (this article) " VLA observations (this article) '^ 



